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FOREWORD 


The  study  reported  here  introduces  a  novel  methodology  to  predict  the 
(frequency-dependent)  effective  material  parameters  characterizing  the  dynamic 
behavior  of  viscoelastic  substances  containing  many  randomly  located  air-filled 
perforations.  These  composite  materials  have  uses  as  underwater  sound-absorbers. 
The  methodology  described  here  is  an  extension  of  earlier  work  of  the  authors 
which  pertains  to  the  case  of  gas-filled  perforations  in  non-absorbing  matrices. 
That  prior  work  was  extended  here  to  the  case  of  absorbing  matrices  containing 
ensembles  of  cavities  of  various  sizes  following  arbitrary  size-distribution 
functions.  The  method  accounts  for  the  effect  of  resonances,  for  arbitrary 
levels  of  viscosity,  for  arbitrary  cavity-size  distributions,  and  it  is  fundamental 
insorfar  as  it  generates  direct  predictions  accounting  for  all  these  effects 
starting  straight  from  the  basic  principles  of  Continuum  Mechanics.  Computer 
codes  to  implement  the  model-predictions  were  generated,  and  a  large  number  of 
pertinent  plots  of  the  frequency  dependence  (at  fixed  concentrations)  or  of  the 
concentration  dependence  (at  fixed  frequencies)  of  the  various  effective  moduli 
and  other  material  descriptors  have  been  computed  and  displayed  in  88  graphs. 

Under  various  conditions  the  present  results  reduce  to  many  of  the  earlier  results 
available  in  the  literature  which  serve  as  check  points.  We  have  identified  four 
relevant  situations  where  additional  work  is  required  to  develop  corrections  to 
this  general  approach. 
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INTRODUCTION 


We  have  developed  a  general  theory  to  model  and  predict  the  effective 
dynamic  material  properties  (i.e.,  moduli,  wavespeeds,  attenuations,  densi¬ 
ties,...)  describing  the  behavior  of  elastic  matrices  containing  .random  distribu¬ 
tions  of  fluid-filled  spherical  cavities.^-  We  have  taken  that  theory, *  which 
accounts  for  resonances  and  frequency-dependent  effects,  and  extended  it  here  to 
the  case  of  lossy  or  sound-absorbing  viscoelastic  matrices  containing  the  same 
types  of  distributions  of  fluid-filled  cavities.  The  methodology  followed  to 
generate  this  extension  is  based  on  scattering  theory,  and  is  the  same  we  used  to 
treat  the  case  of  purely  elastic  matrices.  The  Kelvin-Voigt  viscoelastic  model 
was  used  to  describe  the  rubber  behavior  in  the  manner  we  have  explained 
earlier. 2  A  summary  of  the  basic  results  is  given  below. 


1G.  Gaunaurd  and  H.  ijberall,  "Resonance  Theory  of  the  Effective  Properties  of 
Perforated  Solids,"  J.  of  the  Acoustical  Soc.  America.  Vol.  71,  No.  2,  1982, 
pp.  282-295. 

2 

G.  Gaunaurd,  "Methods  For  Solving  the  Viscoelastic  Equations  For  Cylinder  and 
Sphere  Problems,"  NSWC  TR-76-20,  22  Mar  1976,  28  pp.,  (ADA025302) . 
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I.  THEORETICAL  PREDICTIONS  FOR  LOSSY  MATRICES 


The  effective  properties  are  determined  by  this  model  by  relations  of  the 

type 


which  were  derived  in  ref.  1  (i.e.,  Eq.  (3)). 

The  scattering  coefficients  pertinent  to  the  analysis  to  be  developed  here 

are  A  ,  A,  and  A„.  Coefficient  A  takes  the  form 
o  l  i  o 


A 

o 


x.3 

A 

3i 


,2  2  2 

3  -  x  _  -  x  +  x  . 
of  os  oi 


2  .  2  22 
x,  (1  +  ix,  )  -  x  .  -  x  +  x  . 
d^  d^  of  os  oi 


(2) 


where  x  ,,  x  and  x 

of  os  oi’ 


now  for  lossy  matrices  are 
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We  assume  ChaC  Che  matrix  is  lossy,  while  Che  fluid  inclusions  (medium  2)  are 
noc. 

3 

The  definitions  of  8g^  which  we  have  given  before,  can  now  be  used  to 

express  the  wave  numbers  in  the  viscoelastic  matrix  as  follows. 


ft  +  46.2  +  1  /  /l  +  4 6,2  -  1 

dl  .  ,  2 rf  /l  V  dl 


1  +  48,2 


1  +  48,2 


1  +  48  2  +  1 


1  +  46  2 


Sl  S/1  -  12Ss,  \ 


1  +  48  2  *  c  A/  2  1+46 

S1  S1 V  31 


In  the  absence  of  losses,  these  wavenumbers  reduce  to  the  known  expressions  u>/c 
and  w/cg^,  respectively.  Substituting  Aq  into  the  corresponding  relation 


A  (R)  -  £  A  (a  ) 

0  j-1  °  j 


3g.  Gaunaurd  and  H.  Uberall,  "Theory  of  Resonant  Scattering  From  Spherical 
Cavities  in  Elastic  and  Viscoelastic  Media,"  J.  Acoustical  Society  of  America. 
Vol.  63,  1978,  pp.  1699-1712. 
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yields  Che  following  expression: 


1  - 


plCd2 

"  21sh  > 
p2*d2  dl 


°lcs2 
4  A  S1 

"7  ~  -y  2 ^  “  2i6  ^ 

3  p2*d2  S1 


P,c  2 
4  2  s2 


p; c  2  p,c  2 

4  3  sl  A  2  S1 

1+" - HI  -  2i6 

3  p.c  ,2'  s.  3  P-c ,2 

2  d2  1  2  d2 


=  V 


(6) 


where 


V  = 


i  -  i 


l  - 


P,c,2  0,c  2 

3  dl  4  3  S1 

7Ti(1  -  2iSd  >  -  3  7T2(1  -  2i6s  > 
P2cd.  1  P  P1  S1 


PlCd2 


plCs2 


(7) 


1  -  k  —-4(1  “  2l8a  >x^2(1  +  ix,  )  +  4  —4(1  -  2i6  ) 

3  p2cd2  dl  d  d  3  o2cd2  Sl 


defines  the  parameter  V  and  x^..  =  •'di  ai  •  Manipulating  this  expression  in  the 

/entually  i 


manner  outlined  before  we  eve 


Ly  find 


G  +  iH 


P(a,f)  +  iQ(a,f) 


a  g(a)da 


(8) 


where 


P.c  2  p  c  2 
1  dl  4  1  S1 

1_  P2‘d22  +  3V^  * 


plcd2 


26, 


p  i c  j  2 

4  1  dl 

3  P2cd2 


26s. 


(9a) 
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plCs2 
4  S1 

P(a,f)  =  1  +  ~ 


1^1  (  2rfa  ^  2 


3  P,c,2  3  p„c,2 

2  d2  2  d2 


Q(a,f) 


plcs2 


(9b) 


and  g(a)  is  the  presumably  known  distribution  of  cavity  sizes.  The  real  and 
imaginary  parts  of  4*  are  defined  as 


Vf)  5  T  jf 
3  !f  f 

Jo 


H q(,,» .  a3g(a)da_ 

i  “  I  rs\/  1  ^ 


CP(a,f)]  +  [Q(a,f)r 


HP<a,f)  -  GQ(a,f) 


o  tP(a,f)r  +  [0(a,f )] 


~2  a  g(a)da, 


where  G,  H,  P,  and  Q  are  as  before  and  where 


*F  ■  +  i  '?1  . 


(10a) 


(10b) 


(ID 


We  note  that  4*  (c/o  Eq.  (7))  is  frequency-dependent. 
x<h  =  0),  then  the  ratio  multiplying  (aj^/R3)  in  Eq. 
ana  could  come  out  of  the  sum.  In  this  static  case, 
to  the  concentration  defined  to  be 


If  it  were  not  (i.e.,  if 
(7),  would  not  depend  on  j 
41  would  just  be  proportional 


4> 


(12) 


In  the  dynamic  case  analyzed  here,  this  is  not  the  case  and  the  concentration  is 
defined  to  be 
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4ir 

3  Jo 


g(a)a  da  . 


(13) 


If  all  the  cavities  are  of  the  same  size  a0,  then  all  the  aj  will  be  a0  and  then 


8(a)-  4^3  5(3  "  V* 

TR 


and  either  Eq.  (12)  or  (13)  will  yield  $  =  Na  ^/R^. 
distribution  case  we  have: 


*  -  (?-±jA 

y  +  iq  y 


In  this  uniform  cavity-size 


(14) 


and 


Yf) 


(GP  +  HQ)$ 


[P(aQ,f)]2  +  [Q(ao,f)]2 


(15a) 


^(f)  = 


(HP  -  GQ)  $ 


[P (a  ,f)) 2  +  [Q(a  ,f)]2 
o  o 


(15b) 


So,  we  will  use  either  Eqs.  (10)  or  Eqs.  (15)  depending  on  whether  there  is  a  non- 
uniform  distribution  of  cavity  sizes  or  all  the  cavities  are  of  the  same  sizes, 
respectively. 

w 

The  effective  bulk  modulus  ke  is  now  given  by  the  (normalized)  expression 


plcd 


2 

1 


4 

3 


2iB  ) 
S1 


(16) 


which  can  be  re-written  in  the  following  form 

AS 

k 

e2  a  0_+jLg  m  xu  +  rO  /xv  -  yu\  ,  (17) 

PlCd2  x+iY  x2+y2  Vx2  +  y2/ 
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where  the  parameters  X,  Y,  U  and  V  are  given  by 


X  -  1  -  »R(f)  , 


Y  =  -*T(f) 


(18) 


U  -  1  -  ^ 

3  p 


PlCs2  r 


4 

V  -  T 
3  p 


o.c  2  r 

i  si  r 

picd2  L  Si 


(i  -  v 


28. 


(19) 


(20) 


and  are  given  by  either  Eqs.  (10)  or  (15).  Another  alternative  form  still, 

is, 

(21) 


plCd2 


(1  -  V  +  2  Vi  4  p1cs2 

.(1  -  \)2  +  'i'j2  *  PlCd2 


+  i 


3,  (1  -  ’f  )  p.c  2 

d1  R  Q  Is 

— -  +  #  3  1 


(1  -  *R)2  +  fx2 


3  slPlCd2 


The  expressions  for  the  effective  shear  modulus  and  the  effective  shear  speed 
emerge  from  coefficient  A_.  We  had  already  found  it^  for  elastic  matrices  to  be 


4x,2  p  c  2(p  -  P. c  2) 

1  s^  i  i  s^ 

3i  '  u,(6p  c,2  +  4p,c  2)  +■  pc  2(90,0,2  -  4p  c  2) 

2  1  d^  1  Si  1  si  1  d^  1  si 


(22) 


For  a  viscoelastic  matrix  with  fluid  inclusions  we  find,  in  an  analogous  way, 


4x  3  P.c  2  (1  -  2i8  ) 

dl  1  S1  S1 

3i 


p2  -  PlCg2(l  -  2i8s 


A 


R  +  S 


where 


R  * 


[‘Vdjd  -  12\>  +  “V,2'1  - 
.  >  K=d2 

i  L  i 


S  -  p.c  2(1  -  2i8 
1  s^  s. 


2(1  -  218.  )  -  4p,c  2(1  -  2i8 
d  Is,  s 

1  1 


A 


(23) 


(24a) 


(24b) 
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Substituting  Eq.  (23)  into  the  basic  "effective  property"  relation  for  A 


VR)  "  2  A2(aJ) 
j  -  1 

eventually  yields  the  effective  (normalized)  shear  modulus 


^2  M  +  iN  MK  +  NL  / NK  -  Ml\  ,  (26) 

“1  ‘  K  +  1L  ♦  l*  U2  ♦  L/  . 

where  M,  N,  K  and  L  are 

M  -  f (9p  c  2  -  4p  c  2)  +  26.  (8p,c2  6.  -  18p,c2  6 .  )  1  (1  -  ») 

L  1  “l  1  S1  *1  1  S1  S1  1  “l  dl  J 

H  •  [  <8»1%18»1  -  -  28.1<’»l'dJ  -  '“V,2)]  (1  -  4) 

K  -  (9p  c.2  -  4p  c2  )  +  (6p,c,2  +  4p.c  2)<t> 

1  di  1  Si  1  di  1  Si 

L  -  (8p.c2  &  -  18Plc2  6.  )  -  28„  (9p.c2  -  4p.c*  )<t>  .  (27) 

X  1  a.  X 

Once  the  effective  bulk  modulus  and  the  effective  shear  modulus  are  known,  we  can 
construct  the  effective  dilatational  (or  compressional)  modulus  by  the  relation. 


.  =  p_c  ,2  «.  k  +  -ry_  . 

a  2  ™  ^2  e2  **  * 


Substituting  Eqs.  (26)  and  (17)  into  (28)  eventually  yields. 


2  U  +  iV  XU  +  YV  ( XV  -  YU  \ 

tcd2  ■  X  +  iV  ■  X2  +  ,2  \X2  +  V2j 


where  U  and  V  (without  carets)  are  now 
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a,c  2 

1  -  —  — — — — i 

3»1^ 


1  - 


MK  +  NL 

2  2 
K  +  L 


(1  -  V  )  -  (  28  +  ~  \  T 

R  1  S1  K2  +  L2  1  1 


(30a) 


PlCs2 

V  ,  A— 

3plCd2 


1  - 


MK  +  NL 

2  2 
K  +  L 


(V  +  (1- V 


-  28 


V 


(30b) 


where  X,  Y,  M,  N,  K  and  L  are  as  above.  An  alternative  form  for  this  (normalized) 
effective  dilatational  modulus  is, 


/V  r*  0 

P2Cd3 

1  - 

VR  +  2?^ 

4 

plCs2 

(  MK  +  NL  ^ 

PlCd2  " 

(1  - 

2  2 
*R)  + 

~  3 

p!Cd2 

( 

+  i 


v'V'-v  A  pics? 
1  + 4  1 


(1  -Yr)2  +  ¥x2 


3  PlCd2 


28 

S1  K2  +  L2 


(31) 


These  expressions  for  the  effective  dilatational  modulus  P2cd2  can  s0^ve<^  f°r 
the  effective  dilatational  wavespeed  c^-  The  effective  wavenumber  <d2  *s 


(32) 


and  it  is  therefore  possible  to  relate  the  (real)  effective  sound  speed  C2  and 
the  real  attenuation  <*2 ,  to  the  (complex)  effective  dilatational  speed  c^*  just 
determined,  by  the  following  relations: 
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where  X,  Y,  U  and  V  are  all  given  above.  We  note  that  all  formulas  given  here 
reduce  properly  to  those  given  in  Ref.  1  for  fluid-fillers  inside  the  cavities 
contained  within  elastic  matrices.  Furthermore,  if  the  matrix  is  made  to  be  a 
fluid  matrix  (i.e.,  gas  bubbles  in  a  liquid),  then  all  the  results  here  and  in 
Ref.  1  reduce  properly  to  our  earlier  results  for  bubbly  liquids.^  This  com¬ 
pletes  this  summary  of  theoretical  predictions  for  absorbing  matrices. 

A/  V 

The  effective  Lamd  parameter  X2  can  be  computed  from  the  effective  bulk  k^, 
dilatatlonal  $[^2  ^  \  +  ^~2’  anc*  s^ear  ~2  moc*u-^  hy  the  relations 


4G.  Gaunaurd  and  H.  Uberall,  "Resonance  Theory  of  Bubbly  Liquids,"  J.  of  the 
Acoustical  Soc.  America,  Vol.  §9_,  1981,  pp.  362-370. 
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II.  NUMERICAL  RESULTS 


We  will  concentrate  on  air-filled  (spherical)  cavities  in  a  viscoelastic 
(rubber)  matrix.  We  first  will  consider  the  rubber  and  the  air  material  para¬ 
meters  as  we  have  given  them  in  our  previous  Ref.  (1).  The  rubber  density  is 
”  1.13  g/cm3,  and  that  of  air  is  02  -  0.0012  g/cm^.  The  other  parameters  are 
given  in  Section  III  of  Ref.  (1). 

For  the  viscous  properties  of  the  rubber  we  will  set  8S^  fixed  at  5%  and 

consider  two  values  of  8^,  viz.,  1%  and  5%.  These  quantities  controlling  the 

absorptive  parts  of  and  Ks^  (see  Eqs.  (4))  were  defined  in  our  Ref.  (3).  In 

our  earlier  notation, ^  to  have  8S^  *  5%  means  that  =  10?.  To  have  =  1%  or 

5%  means  that  F  -  4.276  x  10^,  or  F  =  2.21S  x  10^,  respectively.  These  numbers 
were  arbitrarily  chosen  to  introduce  some  level  of  viscosity  in  the  rubber,  and 
to  see  how  the  non-absorbing  results  plotted  in  Ref.  (1)  were  altered  by  the 
presence  of  matrix  viscosity.  However,  these  values  (8d1/''  5%,  8Sl~  5%)  are  not 

uncommon  in  practice.  We  have  chosen  concentrations  of  <t>  =  5  and  10%  because 
these  were  the  numbers  selected  in  Ref.  (1),  and  it  was  of  interest  to  see  how  the 
"effective  properties"  at  these  concentrations  were  altered  from  the  values  in 
Ref.  (1),  to  the  present  values  in  the  presence  of  matrix  losses.  We  have  stated 
in  earlier  work  of  ours  that  the  present  theory  neglects  re-scattering  (i.e., 
multiple  scattering)  from  the  many  cavities  in  the  matrix.  This  assumption 
ceases  to  hold  for  concentrations  above  25%,  which  is  the  upper  limit  of  validity 
of  this  theory,  without  additional  multiple-scattering  corrections.  In  the  first 
few  plots  that  follow,  the  cavities  are  of  uniform  sizes  as  we  will  describe  next. 

Figure  1  shows  the  real  and  imaginary  parts  of  the  (normalized)  effective 
dilatational  modulus  Mj,  versus  frequency  for  clouds  of  cavities  of  2um  in  radius 
at  5%  and  10%  concentrations  when  the  matrix  has  a  viscous  level  of  Bd^  *  5%,  and 
Bs^  *  5%.  The  formula  plotted  is  Eq.  (29).  For  the  materials  and  concentrations 
considered,  these  plots^are  undistinguishable  from  the  corresponding  plots  for  the 
effective  bulk  modulus  ke^  as  given  by  Eq.  (17). 


Figure  2  shows  the  same  quantities  at  the  same  concentrations  and  viscous 
levels  but  for  clouds  of  larger  cavities  all  at  uniform  20  micron  sizes.  The 
resonance  effects  are  quite  visible  and  are  clearly  shifted  to  the  lower  fre¬ 
quency  region.  The  formulas  plotted  are  Eqs.  (29)  and/or  (17),  which  yield  un¬ 
distinguishable  results  for  the  present  materials.  At  the  high  frequency  ends 
of  the  real  plots  (left)  we  see  that  Md2  •+■  Md^,  which  means  that  the  effective 
properties  coincide  with  those  of  the  matrix  at  high-frequencies  and  low 
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I 


concentrations.  The  numbers  in  the  abscisas  are  frequencies  in  megahertz. 

The  results  in  Figs.  1  and  2  are  to  be  compared  to  those  in  ref.  (1),  which 
were  computed  for  the  same  parameters  but  without  absorption  in  the  matrix  (i.e., 
6^  *  0,  Ss^  “  0) .  Viscosity  shifts  the  location  of  the  (giant  monopole) 
resonance  to  lower  frequencies  and  knocks  down  its  amplitude  in  the  real  plots. 

In  the  imaginary  plots,  the  amplitude  of  the  dip  is  less  pronounced  the  higher 
the  viscosity,  and  the  location  of  the  dip  is  also  shifted  to  the  low  frequency 
end.  Viscosity  also_broadens  all  these  resonance  peaks  (or  dips)  so  that  the 
widths  are  all  increased  with  absorption,  as  it  should  be.. 

Figures  3  and  4  are  companion  figures  to  each  other  and  they  display  the 
(normalized)  effective  sound  speed  and  attenuation  for  a  cloud  of  air-filled  ca¬ 
vities  in  lossy  rubber  at  $  »  5%  (upper  graphs)  and  10%  concentration  (lower 
graphs).  All  the  cavities  are  of  the  same  uniform  size,  aQ  =  2  microns.  Figure  3 
is  for  a  viscous  level  (i.e.,  Bd,  *  1%,  8Sl  *  5%)  and  Fig.  4  is  for  a  higher  level 
(i.e.,  Bd-,  =  5%,  Ss.  *  5%).  Both  these  figures  can  be  compared  to  analogous  ones 
we  have  published  elsewhere^  for  totally  non-absorbing  matrices  (i.e.,  Bd.  *  Bs-.  = 
0) .  Higher  absorptions  shift  the  resonances  toward  the  lower  frequency  end  of  tih.e 
spectrum  and  broaden  the  peaks.  We  have  examined  these  effects  in  the  past^  in 
analogous  instances  for  a  single  cavity.  These  effects  continue  to  be  present  for 
clouds  of  cavities  as  we  can  see.  The  formulas  being  plotted  in  Figs.  3  and  4  are 
given  in  Eqs.  (33)  and  (34).  High  and  low  frequency  limits  could  be  easily  found 
from  those  expressions,  but  these  will  not  be  investigated  analytically  here. 
Numerical  values  are  observable  in  the  plots  for  these  limits. 

Figures  5  and  6  are  companion  figures  to  each  other  and  they  display  the 
(normalized)  effective  sound  speed  and  attenuation  versus  frequency  for  a  cloud 
of  air-filled  cavities  in  lossy  rubber  now  all  of  larger  but  uniform  sizes  (a  = 

20  microns).  The  upper  graphs  are  for  $  =  5%  and  the  lower  ones  for  $  =  10%. 

Fig.  5  is  for  a  level  of  viscosity  (i.e.,  Bd^  =  1%,  BSl  =  5%)  and  Fig.  6  is  for 
a  higher  level  (i.e.,  8d^  ”  5%  ■  Bg, ) .  Both  these  figures  can  be  compared  to 
similar  ones  for  no  viscosity  at  all  in  the  matrix  which  were  published  else¬ 
where.  The  pair  of  Figs.  3  and  4  are  analogous  to  the  pair  formed  by  Figs. 

5  and  6,  the  difference  being  that  the  first  pair  corresponds  to  cavities  2ym 
in  size  while  the  second  is  for  20  micron  cavities. 

We  recall  that  the  effect  of  the  giant  monopole  resonance  for  a  single  gas- 
filled  cavity  in  rubber?  manifests  itself  at  k<j,  Sq  5  0.1,  where  kd^  =  2irf/cd-^. 

For  rubber  of  c^^  ■  1.4  x  10^  cm/sec,  that  condition  translates  to 

fa0  *  2.5  (38) 

^G.  Gaunaurd  and  H.  Uberall,  "Errata:  Resonance  Theory  of  the  Effective  Proper¬ 
ties  of  Perforated  Solids,"  J.  of  the  Acoustical  Soc.  America  (to  be  published), 
1983. 

**G.  Gaunaurd,  et  al.,  "New  Method  To  Determine  Shear  Absorption  Using  the  Visco- 
elastodynamic  Resonance  Scattering  Formalism,"  J.  of  the  Acoustical  Soc.  America. 
Vol.  64,  1978,  pp.  1211-1212. 

^G.  Gaunaurd,  et  al. ,  "Giant  Monopole  Resonances  in  the  Scattering  of  Waves  From 
Gas-filled. Spherical  Cavities  and  Bubbles,"  J.  of  the  Acoustical  Soc.  America, 
Vol.  65,  1979,  pp.  573-594.  - L 
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where aQ  is  in  centimeters  and  f  is  in  kilohertz.  When  there  are  many  cavities 
in  the  rubber,  the  number  on  the  right  side  of  Eq.  (38)  is  a  bit  higher  (c/o 
Figs.  3-6)  but  we  will  use  Eq.  (38)  for  the  argument  that  follows.  If  we  have 
f  ■  20  kHz,  the  first  resonance  occurs  for  cavities  of  size  a0  =  0.125  cm.  If 
f  =  60  kHz,  then  the  required  cavity  size  for  resonance  is  ao  =  0.0417  cm.  So, 
cavity  sizes  required  for  resonance  in  this  frequency  range  go  from  0.4  mm  to 
1.2  mm. 


The  effective  shear  modulus  is  displayed  versus  concentration  <t> 

(0  <  $  <  50%)  at  three  levels  of  viscosity  in  Fig.  7.  These  graphs  exhibit  the 
real  and  the  imaginary  parts  of  the  normalized  modulus,  all  computed  for  air- 
filled  cavities  in  rubber.  The  equation  being  plotted  is  Eq.  (26).  For  the 
three  absorption  levels  shown  in  Fig.  7,  the  real  parts  of  the  modulus  exhibit 
very  small  changes  while  the  changes  in  the  imaginary  parts  are  more  perceptible. 

Figure  8  shows  the  real  and  imaginary  parts  of  the  (normalized)  effective 
shear  speed  versus  concentration  $  for  two  levels  of  viscosity  (i.e.,  3<j^  =  1%, 

Ss^  =  5%  and  8^  =  8S^  =  5%  as  they  are  obtained  from  Eq.  (36).  Higher  absorp¬ 
tions  reduce  the  shear  speed  in  the  perforated  solid  and  also  its  corresponding 
imaginary  part. 

A  relation  for  the  effective  density  such  as  Eq.  (35)  was  plotted,  in 
normalized  form  in  earlier  work  of  ours  (viz.,  ref.  (1),  Fig.  6)  and  will  not 
be  repeated  here. 

The  computer  programs  developed  to  generate  Figs.  1-8  could  also  be  used  to 
plot  any  of  the  formulas  in  Section  I,  for  any  material  combination,  absorption 
levels,  concentrations,  or  distributions  of  cavity  sizes  and  frequency  ranges. 
Appendix  1  shows  the  complete  listings  for  these  programs. 

These  programs  have  been  used  to  compute  Figs.  9-12.  Figure  9  shows  the 
real  (left)  and  imaginary  (right)  parts  of  the  (normalized)  effective  bulk 
modulus  versus  frequency  for  a  uniform  distribution  of  cavity  sizes  all  of 
radius  aQ  =  0.125  cm  at  <t  =  10%.  The  upper  plots  are  for  a  viscous  level  of 
Bd-j^  *  Bs^  *  5%  while  the  lower  plots  are  for  non-viscous  matrices  (Bd.  =  Bg.  =  0^ 
in  order  to  show  the  effect  of  (light)  matrix  absorption.  As  in  previous  figures, 
these  results  are  numerically  equal  to  those  of  the  (normalized)  effective  dila- 
tational  modulus.  The  pertinent  equations  plotted  are  Eqs.  (17)  or  (29).  The 
frequency  ranges  are  now  restricted  to  the  band  0  <  f  <  40  kHz. 

Figure  9  shows  the  effective  (and  normalized)  sound  speed  and  attenuation 
versus  frequency  plots  for  clouds  of  cavities  of  (uniform)  size  a0  =  0.125  cm 
at  <t>  «  10%.  The  upper  plots  are  to  be  compared  ot  the  lower  plots,  to  see  the 
effect  of  viscosity.  The  upper  graphs  are  for  the  viscous  level  Bdi  ■  Bs,  "  5%, 
while  the  lower  ones  are  for  non-lossy  matrices  (i.e.,  Bd^  “  6Sl  -  6%)  and  are 
shown  for  comparison.  Resonance  effects  are  clearly  observable  around  20-24 
kHz  as  one  would  expect  from  Eq.  (38)  for  this  size  cavities,  which  are  the  ones 
resonating  for  incident  waveforms  of  frequencies  20  kHz.  Observation  of  these 
plots  show  that  viscosity  in  the  matrix  broadens  the  resonance  peaks  and  shifts 
them  to  higher  frequencies.  A  more  detailed  analysis  of  this  effect  could  be 
done  with  results  such  as  these.  The  formulas  plotted  in  Fig.  10  are  Eqs.  (33) 
and  (34). 
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Figures  11  and  12  are  the  identical  counterparts  of  Figs.  9  and  10, 
respectively,  but  for  clouds  of  cavities  of  size  a0  =  0.0417  cm.  These  are  the 
sizes  that  excite  the  giant  monopole  resonance  (at  around  70  kHz)  whenever  the 
incident  wave  is  of  frequency  f  =  60  kHz  (c/0  Eq.  (38)).  The  frequency  range  in 
these  figures  now  includes  the  broader  band  0  ^  f  ^  100  kHz.  The  broadening  and 
shift  of  the  attenuation  (i.e.,  “a 2 )  peak  with  increased  absorption  is  clearly 
visible  in  Fig.  12.  All  plots  presented  in  Figs.  1-12  are  for  clouds  of  air- 
filled  cavities  in  rubber  having  uniform  sizes.  In  these  cases,  and  are 
always  found  from  Eqs.  (15). 

It  is  of  interest  to  see  how  these  quantities  and  fj,  as  computed  from 
Eqs.  (15),  behave  as  functions  of  frequency.  Figure  13  displays  these  quantities 
in  the  range  0  ^  f  ^  40  kHz  for  swarms  of  air-filled  cavities  of  uniform  size  a0  * 
0.125  cm,  and  a  $  ■  10%  concentration.  The  upper  graphs  are  for  no  matrix 
viscosity,  while  the  lower  ones  are  for  a  moderate  matrix  viscosity  of  Bd^  *  Bs^  * 
5%.  It  is  quite  clear  that  viscosity  shifts  and  broadens  the  resonance  peaks. 

The  rest  of  the  figures  (i.e..  Figs.  14-22)  are  all  for  the  case  where  the 
cavity  sizes  are  not  uniformly  distributed.  The  methodology  presented  here  can 
analyze  and  account  for  arbitrary  size  distributions  g(a).  We  have  selected  a 
triangular  distribution  of  cavity  sizes  to  illustrate  this  effect  in  a  simple 
case,  although  real  cavity  size-distributions  are  more  likely  to  be  Gaussian, 
provided  the  cavities  are  really  randomly  located,  as  opposed  to  periodically 
located  in  a  lattice.  Figure  14  shows  the  distribution  function  and  the  cummula- 
tive  distribution  function  of  the  triangular  distribution.  The  peak  of  the  dis¬ 
tribution  is  set  to  occur  at  value  a0.  It  can  be  easily  verified  that  the  mean 
of  the  triangular  distribution  is  1.5  times  the  mean  of  the  uniform  distribution 
function. 

Figure  15  shows  the  same  quantities  plotted  in  Fig.  13  (viz.,  ^r  and  fj)  in 
the  absence  (top)  and  presence  (bottom)  of  viscosity  in  the  matrix,  at  a  10% 
concentration  level.  The  cavities  are  now  triangularly  distributed  with  the  peak 
of  the  triangular  distribution  a0  at  0.125  cm.  Comparison  of  Figs.  13  and  15 
shows  that  for  triangularly  distributed  cavities,  the  resonance  peaks  are  shifted 
to  lower  frequencies  and  broadened  to  wider  widths.  The  magnitude  of  the  peaks 
is  also  reduced.  So,  cavity  size-distributions  other  than  uniform  have  an  effect 
similar  to  that  of  viscosity  in  the  matrix,  viz.,  they  lower  and  broaden  the 
resonances,  and  knock  down  their  respective  amplitudes.  Figure  16  is  identical 
to  Fig.  15  except  that  it  has  now  been  recomputed  for  the  higher  concentration 
level  of  »  15%.  These  observations  all  emerge  from  quick  inspections  of  Figs. 
13  and  15,  and  they  could  be  quanitatively  analyzed  in  further  detail.  Figure  16 
shows  that  a  higher  concentration  of  triangularly  distributed  cavities  further 
increases  the  size  of  the  peaks,  all  other  things  being  equal. 

Figure  17  is  the  counterpart  of  Fig.  9  for  triangularly  distributed  swarms 
of  air-filled  cavities  in  a  rubber  matrix  at  10%  concentrations.  Now  a0  -  0.125 
cm  is  the  peak  of  the  triangular  distribution.  The  upper  graphs  are  for  moderate 
matrix  viscosity  (i.e.,  $d^  *  Bg,  *  5%)  while  the  lower  ones  are  for  no  viscosity 
at  all.  This  Fig.  17  displays  the  real  and  the  imaginary  parts  of  the  effective 
bulk  modulus,  versus  frequency  in  the  range  0^f  ^40  kHz.  The  graphs  in  Fig.  17 
are  very  close  to  those  in  Fig.  9,  except  in  the  middle  portion  (f  **  16  kHz) 
where  an  additional  dip  is  observable.  These  are  the  regions  where  Yr  and  fj 
showed  a  difference  (Fig.  15)  in  relation  to  their  values  for  uniform  cavity-size 
distributions  (c/o  Fig.  13) 
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Figure  18  is  the  counterpart  of  Fig.  10,  now  for  triangularly  distributed 
clouds  of  cavities  at  10%  concentrations,  where  a0  »  0.125  cm  is  the  peak  of  the 
triangular  distribution.  Figure  18  displays  in  this  case  the  effective  sound 
speed  and  effective  attenuation  a2a0  versus  frequency  plots.  Comparison 

of  Figs.  18  and  10  shows  that  the  triangular  distribution  of  cavities  has  a  very 
marked  effect  on  the  attenuations.  The  effective  sound  speed  also  changes — it  is 
lower  now  at  all  frequencies — and  there  is  a  smoothing  effect  in  the  mid-portions 
of  the  graphs  where  the  resonances  in  changed  substantially.  The  plots 

are  computed  at  two  levels  of  viscosity  (viz.,  =  Bs^  =  5%,  upper  graphs,  and 
*  8gl  =*  0%,  lower  graphs)  as  before. 

The  results  in  Fig.  17  have  been  re-computed  in  Fig.  19  for  a  15%  concentra¬ 
tion,  all  else  being  the  same.  Higher  concentrations  tend  to  flatten  these 
spectra,  and  we  indeed  see  that  at  low  frequencies  the  real  and  imaginary  parts 
of  the  effective  bulk  moduli  in  Fig.  19,  computed  assuming  a  triangular  cavitj - 
size  distribution,  are  lower  than  in  Fig.  17  where  the  concentration  was  only 
10%.  At  the  higher  frequency  end  of  these  spectra,  we  observe  higher  values  in 
Fig.  19  than  those  in  the  corresponding  plots  of  Fig.  17.  Thus,  higher  concentra¬ 
tions  have  clearly  flattened  the  earlier  results. 

Figure  20  is  the  same  as  Fig.  18  but  re-computed  at  15%  concentrations.  It 
is  clear,  from  inspection  of  these  figures,  that  higher  concentration  levels  tend 
to  lower  the  effective  sound  speeds  at  all  frequencies  and  to  increase  the 
effective  attenuations,  particularly  at  the  high  frequencies  beyond  that  of  the 
giant  monopole  resonance  (^20  kHz  in  this  case). 

We  have  recomputed  Fig.  19  for  smaller  cavities,  still  triangularly 
distributed,  but  now  with  peak  at  a0  ”  0.0417  cm.  The  result  is  shown  in  Fig.  21. 
This  figure  shows  the  real  and  the  imaginary  parts  of  the  effective  bulk  modulus 
versus  frequency  (0  <  f  <  100  kHz)  plots,  for  ensembles  of  triangularly  distri¬ 
buted  cavities  in  rubber,  at  15%  concentrations.  The  upper  graphs  are  for  a 
viscous  level  Bd-^  *  Bs-^  *  5%,  while  the  lower  graphs  are  for  Bd^  “  8Sl  =  0%,  (no 
viscosity).  Figure  21  (triangular  distribution)  can  also  be  compared  to  Fig.  11 
(uniform-size  distribution),  although  the  concentrations  are  not  the  same.  The 
same  general  remarks  made  before  still  hold  in  this  case. 

Figure  22  is  the  re-computation  of  Fig.  20  but  for  smaller  a0.  Therefore, 
Fig.  22  gives  the  effective  sound  speed  and  attenuation  versus  frequency  (0  f 
<  100  kHz)  plotts,  for  ensembles  of  gas-filled  cavities  in  a  rubber  matrix  all 
triangularly  distributed  with  the  peak  a0  now  at  0.0417  cm.  The  concentration 
level  is  15%  and  the  same  two  earlier  viscous  levels  (Bj.  *  Bs.  *  5%  and  Bd-^  * 

Bs^  *  0%)  are  considered.  Figure  22  can  also  be  compared  to  Fig.  12  although  the 
concentrations  are  slightly  different  (i.e.,  10%  versus  15%).  Again,  the 
triangular  distribution  lowers  the  effective  sound  speed,  particularly  at  high 
frequencies,  and  increases  the  effective  attenuation  over  a  wider  frequency  band, 
mainly  because  of  the  similar  flattening  and  broadening  effects  non-uniform 
distributions  and  higher  viscosity  levels  have  on  the  resonance  peaks.  All  these 
effects  are  quantitatively  exhibited  in  these  22  figures  which  contain  88  graphs, 
all  computer  generated,  from  the  basic  theory  presented  here. 
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III.  DISCUSSION  OF  THE  CHABAN  RESULTS 


The  Chaban  results  are  valid  for  non-absorptive  matrices  and  they  are 
frequency-independent  (i.e.,  they  are  static  results).  Their  derivation  is  based 
in  the  so-called  self-consistent  approach  which  has  been  the  subject  of  much 
controversy  and  doubt. 


For  non-absorptive  rubber  matrices  (i.e.,  =  0,  6S^  =  0)  containing 

clouds  of  uniform  size  air-filled  cavities  of  radii  aQ,  the  (normalized)  effective 
bulk  modulus  given  in  our  Eq.  (17)  reduces  to  the  expressions  given  in  Eqs.  (34)- 
(40)  of  ref.  (1).  In  the  static  or  long-wavelength  limit  (i.e.,  f  -*■  0) ,  that 
expression  contains  the  well-known  Kemer  result^ 


ke.  +  Tl 

kia-#)+SlCTiT,# 


ke.  +  Tl 

1  ■  * +l  — ■ ~ 4ir  14 

e2  3  1 


which  for  dilute  concentrations  ($  «  1)  yields  the  expression 


k  -  k 


[(v*0/  N*^)]N'k*») 


which  is  Chaban's  dilute  limit3  (ref.  9,  Eq.  (41). 


(39) 


(40) 


The  (normalized)  effective  shear  modulus  given  in  Eq.  (26)  reduces  to  the 
expression  in  Eq.  (45)  of  our  ref.  1  when  there  are  no  losses  in  the  rubber 
matrix  (8dl  *  Bg.  ■  0).  This  later  equation,  in  the  dilute  limit  ($  «  1)  agrees 
with  Chaban's  dilute  limit  formula  which  is 


Pi  -  |  *(1  -  —)(*!  +  Zy^E, 


(41) 


Q 

E.  H.  Kerner,  "The  Elastic  and  Thermoelastic  Properties  of  Composite  Media," 
Proc.  Phya.  Soc.  (London),  Vol.  69B,  1956,  pp.  808-813. 

I.  A.  Chaban,  Soviet  Phys.  -  Acoustics,  Vol.  H),  1965,  pp.  298-304,  and,  Vol.  11, 
1965,  pp.  81-86. 
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where  E  is  a  parameter  defined  by 


,  U2 


31^  +  8y^  9A^  +  14y^ 


(42) 


It  is  not  hard  to  show  that  this  expression  coincides  with  the  expression  we  gave 
in  our  ref.  1  (Eq.  (80)),  viz.. 


in  terms  ofthe  Poisson's  ratio  the  matrix.  Therefore,  our  frequency- 

dependent  method  presented  here  for  absorbing  matrices,  greatly  extends  and  con¬ 
tains  the  static  and  non-absorbing  results  of  Chaban  in  the  dilute  limit. 

Finally,  for  air-filled  (or  evacuated)  cavities  in  rubber  (assumed  lossless) 
Eqs.  (39)  and  (43)  yield,  approximately. 


(44) 


(45) 


which  are  results  due  to  Gassmann  and  Dewey,  respectively,  and  which  we  already 
studied  in  our  earlier  work.^-  These  static  formuli  are  sometimes  quoted  by 
absorber  designers. 
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IV.  DISCUSSION  OF  ASSUMPTIONS  AND  REGIONS  OF  VALIDITY 
OF  THE  METHODOLOGY 


FUTURE  PLANS 


A  variety  of  assumptions  were  made  during  the  construction  of  the 
methodology  described  here.  Some  of  the  most  important  ones  are  as  follows:  a) 
The  cavities  are  assumed  small  compared  to  the  wavelength  of  the  incident  wave. 

In  spite  of  this  small  a/X  assumption,  resonance  effects  are  present,  and  are 
taken  into  account  by  the  method,  since  the  dominant  giant  monopole  resonance 
occurs  at  very  low  a/X-values.  b)  The  volume  concentrations,  $,  of  cavities 
within  the  matrix  are  assumed  small  (i.e.,  4>  <  25%).  This  implies  that  the 
possible  multiple  scattering  effects  that  would  be  present  for  densely  packed 
configurations,  are  of  secondary  importance  and  to  first  orders,  may  be  neglected, 
c)  All  the  observations  are  made  in  the  far-field  of  the  effective  sphere* 

(i.e.,  r  »  R)  so  that  only  the  global  effects  of  the  ensemble  of  cavities  are 
noticed  at  the  observation  point.  This  condition  is  met  for  either  low  concen¬ 
trations  (i.e/,  $  <<  1),  low  frequencies,  or  small  cavity  sizes  (a/X  <<  1).  d) 
The  giant  monopole  resonance  has  a  dominant  effect  over  all  others  and  manifests 
itself  very  strongly  at  very  low  a/X  -  ratios.  We  have  shown  this  to  be  true? 
for  air-filled  cavities  in  rubber  matrices,  weather  the  rubber  is  assumed  lossy 
or  not.  This  assumption  does  not  hold  true  for  solid  inclusions,  and  the  whole 
analysis  presented  here  must  be  modified  in  that  case,  e)  The  cavities  are  all 
contained  in  a  boundless  matrix  of  infinite  extent  (i.e.,  no  boundaries  are 
present).  Therefore,  the  present  model  works  best  when  the  cavities  are  small 
compared  to  the  layer  thickness.  When  the  incident  waveforms  are  of  low 
frequencies,  the  cavity  sizes  required  to  excite  their  giant  monopole  resonance 
(i.e.,  the  lowest),  could  be  large  enough  to  make  them  comparable  in  size  to  the 
thickness  of  the  matrix  (or  layer)  that  contains  them.  In  this  case,  an 
additional  correction  must  be  introduced  into  the  analysis  to  account  for  the 
nearby  presence  of  boundaries  or  layer  edges.  At  higher  frequencies  of  the  in¬ 
cident  wave,  the  required  cavity  sizes  became  small  compared  to  the  matrix  thick¬ 
ness,  the  thickness  correction  becomes  negligible,  and  the  boundless  matrix 
approach  presented  here  remains  valid.  We  have  described  the  case  where  the 
cavity  size  and  layer  thickness  are  of  comparable  sizes,  as  a  "correction"  to 
the  boundless  matrix  approach  described  here.  This  is  so  because  the  cavity 
size  will  still  control  the  location  of  the  resonance,  however,  the  finite  layer¬ 
thickness  will  shift  it,  broaden  it,  and  modify  it  from  its  unbounded  value.  A 
numerical  example  will  quantitatively  illustrate  these  points. 


Assume  there  is  a  (thick)  rubber  layer  covering  an  underwater  structure  that 
is  10  cm  (4  inches)  thick.  This  is  not  an  unrealistic  thickness  if  the  coating 
is  to  absorb  effectively  in  the  low-end  of  the  frequency  spectrum.  At  fre¬ 
quencies  of  f  =  20  kHz,  we  can  obtain  an  estimate  for  the  cavity  size  required  to 
excite  the  first  (i.e.,  the  giant  monopole)  resonance  from  Eq.  (38).  This 
estimated  value  is  a  =  0.125  cm  (i.e.  one-and-a-quarter  millimeter).  This  size 
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is  about  a  fiftieth  of  the  layer  thickness  and  hence,  the  theory  developed  here 
applies  quite  adequately.  In  fact,  for  such  frequencies,  and  such  (small)  cavity 
sizes,  one  does  not  need  such  a  thick  layer.  In  fact,  a  thickness  of  about  1  cm 
will  still  be  adequate  to  permit  the  application  of  the  present  model  without 
any  additional  corrections.  For  incident  frequencies  of  f  =  60  kHz,  Eq.  (38) 
would  yield  resonating  cavity  sizes  of  value  aQ=  0.042  cm,  which  are  even  smaller 
when  compared  to  the  layer  thickness  than  before,  and  thus  the  layer  thickness 
could  be  reduced  even  further  to  perhaps  1/3  cm.  without  modifications.  If,  on 
the  other  hand  we  had  lower  incident  frequencies  (say,  f  '  2.5  kHz),  then  Eq.  (38) 
would  yield  cavities  of  diameter  2  cm,  which  is  still  relatively  small  compared 
to  the  10  cm  layer-thickness,  provided  there  is  only  one  such  row  of  cavities 
within  the  layer  as  it  is  often  the  case.  It  is  for  these  "large"  cavities 
which  lie  close  to  the  layer  edges  that  we  should  introduce  edge-corrections  or 
f inite-thickness  corrections,  not  accounted  for  in  the  present  analysis. 

There  seems  to  be  no  upper  limit  on  the  levels  of  viscosity  (compressional 
or  shear)  that  could  be  used  within  this  model.  Either  type  of  viscosity  could 
well  be  as  high  as  100%  without  any  foreseeable  complications.  We  have  already 
seen  that  for  single  cavities,  higher  levels  of  viscosity  produce  shifts  in  the 
resonance  frequencies  and  broadenings  of  the  resonance  peaks. ^  This  phenomenon 
has  not  been  examined  in  detail  within  the  context  of  the  present  model  for  many 
cavities  within  a  matrix. 

Inequalities  of  some  design  importance  are  X  >>  L  >>  a,  where  L  is  the  layer 
thickness  and  "a"  the  cavity  radius.  A  further  rule  of  thumb  between  layer  and 
cavity  dimensions  is  L  >  6a. 

The  model  assumes  spherical  cavities.  The  shape  of  the  cavities  appears  to 
be  irrelevant  provided  that  the  air  contents  (i.e,  volume)  of  any  two  possible 
shapes,  as  well  as  their  concentration,  remains  the  same.  This  follows  from 
conditions  a)  and  b) ,  since  we  are  dealing  with  small  cavities  compared  to  the 
layer  thickness  and  the  wavelength  of  the  incident  sound  wave.  Under  the  present 
assumptions  a)  -  e) ,  this  model  will  apply  to  specific  absorber  designs  such  as 
AA-17  provided  one  develops  the  finite-thickness  corrections  for  the  layer  which 
is  needed  under  assumption  e) .  We  make  this  statement  specifying  this  correction, 
even  though  the  frequencies  at  which  AA-17  is  most  effective  are  quite  high, 
because  the  specific  design  calls  for  cavities  which  are  very  close  to  the  layer 
edges.  At  the  relatively  low  (under  25%)  cavity  concentrations  present  in  AA-17, 
it  seems  reasonable  to  assume  that  the  cylindrical  cavities  could  be  replaced 
by  equivalent  spherical  ones  of  the  same  volume  and  concentration,  without  intro¬ 
ducing  much  error  in  the  final  results.  This  conjecture  should  also  be  further 
investigated  and  to  this  end  we  have  already  devised  a  plan  to  study  the  monopole 
resonance  of  finite-length  cylindrical  cavities. 

In  summary,  four  further  corrections  should  be  included  in  the  analysis  for 
various  instances.  For  densely  packed  inclusion  configurations,  additional 
multiple  scattering  corrections  have  to  be  developed.  For  the  low  frequencies 
requiring  the  large  resonating  cavities  having  walls  near  to  the  layer  edges,  we 
must  introduce  corrections  to  account  for  the  finite  layer-thickness.  For 
cavities  with  shapes  other  than  spherical,  particularly  at  high  frequencies, 
shape  corrections  should  be  considered.  Of  particular  interest  is  the  case  of 
finite  cylindrical  cavities,  for  which  even  the  giant  monopole  resonance  is  not 
well  understood.  Fourth,  some  absorber  designs  call  for  the  introduction  of 
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solid  inclusions  in  the  matrix  rather  than  fluid  ones  (i.e.,  gas-filled  cavities). 
This  final  instance  requires  further  study,  for  both  elastic  and  viscoelastic 
matrices.  Finally,  there  seems  to  be  no  upper  viscosity  level  beyond  which  the 
present  model  fails.  We  have  developed  research  plans  to  thoroughly  study  and 
understand  all  four  of  these  modifications  and/or  corrections,  in  future  work. 

The  model  also  takes  into  account  arbitrary  cavity  size  distribution  functions. 
Numerous  calculations  have  been  performed  and  displayed  in  two  basic  instances. 
The  first  is  for  the  case  when  all  the  cavities  are  of  the  same  uniform  size  aD. 
The  second  is  for  the  case  when  some  cavities  are  larger  than  others,  the  whole 
collection  following  the  simple  (but  arbitrarily  chosen)  triangular  distribution 
function  of  peak  (mean)  value  a0.  We  have  produced  many  quantitative  plots  for 
the  variations  of  the  most  significant  quantities,  in  both  instances.  The  over¬ 
all  effect  of  non-uniform  (triangular  in  this  case)  cavity  size  distributions  is 
analogous  to  that  of  a  higher  viscous  absorption  level,  viz.,  it  shifts  the 
actual  resonance  locations,  broadens  their  peaks  and  smoothes-out  or  diminishes 
their  spectral  amplitudes. 
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V.  SUMMARY  AND  CONCLUSIONS 


In  the  present  work  we  have  developed  a  general  theory  to  predict  the  dynamic 
effective  material  properties  of  lossy  perforated  matrices  with  fluid-fillers. 
These  composites  are  considered  to  be  viscoelastic  (i.e.,  sound  absorbing)  matrices 
containing  random  distributions  of  fluid-filled  spherical  cavities  or  perfora¬ 
tions  of  various  sizes.  The  method  accounts  for  the  losses  in  the  rubberlike 
material  in  a  fundamental  way  which  emerges  from  t^e  basic  principles  of 
mechanics  and  the  Kelvin-Voigt  viscoelastic  model.  Hence,  the  method  does  not 
proceed  in  an  ad-hoc  way  as  it  has  been  done  by  authors  who  have  derived^  (or 
repeated  the  derivations  of  earlier-^)  lossless  theories  and  "extended  them"  by 
assuming  complex  wavenumbers  or  complex  wavespeeds  which  can  not  be  related  to 
observed  material  properties  in  a  systematic  (or  any  other)  way. 10  Our  method 
is  based  on  scattering  theoryl  and  thus,  it  yields  frequency-dependent  results 
in  a  natural  way.  This  frequency-dependence  of  results  is  again  found  from  the 
fundamental  laws  of  mechanics  and  it  is  to  be  contrasted  with  other  approaches 
which  are  static  but  then  they  are  "made  to  be  frequency-dependent"  by  intro¬ 
ducing  experimentally  observed  data,  fitted  by  approximate  means,  into  the 
analytic  expressions  for  the  static  results.  Thirdly,  our  method  accounts  for 
the  actual  presence  of  many  cavities  of  many  sizes  within  the  matrix,  the  effects 
of  each  one  added  by  the  technique.  This  is  unlike  other  approaches  that  have 
dealt  only  with  one  such  cavity  in  the  matrix^ » 10  and  merely  stated  that  if 
there  had  been  many,  the  effects  would  have  been  additive  under  some  sort  of 
addition  law  of  effective  behavior  that  (incredibly)  was  never  really  produced. 10 
Although  these  authors^  talk  about  the  effective  behavior  of  the  composite,  the 
composite  they  have  considered  has  only  one  hole,  and  matrix  viscosity  is  intro¬ 
duced  fictitiously  in  the  ad-hoc  fashion  we  mentioned  above. 

Our  method  stays  away  from  the  procedures  used  in  the  "self-consistent 
method"^  because  of  the  conceptual  difficulties  and  the  inconsistencies 


^R.  L.  Kligman,  et  al. ,  "Effective  Dynamic  Properties  of  Composite  Viscoelastic 
Materials,"  J.  of  the  Acoustical  Soc.  America,  Vol.  _70,  1981,  pp.  1437-1444. 

^N.  Yamakawa,  "Scattering  and  Attenuation  of  Elastic  Waves,  Parts  I  and  II," 
Geophys.  Magaz.  (Tokyo),  Vol.  31^,  1962  ,  pp.  63-103. 


‘E.  Meyer,  et  al. ,  "Pulsation  Oscillations  of  Cavities  in  Rubber,"  J.  of  the 
Acoustical  Soc.  America,  Vol.  30,  1958,  pp.  1116-112/'. 
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associated  with  that  approach, 13  that  have  seriously  questioned  its  validity. 

The  present  analysis  yields  results  that  reduce  properly  to  every  previously 
published  static  result^-  not  using  the  self-consistent  method,  and  also  to  the 
few  dynamic  results  that  are  available  for  simpler  matrix/filler  combinations.4 

The  present  method  yields  effective,  frequency-dependent,  analytic  expres¬ 
sions  for  the  effective  moduli  (bulk,  dilatational,  and  shear),  effective  density, 
effective  sound  speed  and  effective  attenuation,  for  a  viscoelastic  (lossy) 
medium  of  large  extent,  containing  many  fluid-filled  spherical  cavities  of  various 
sizes.  All  results  are  checked  by  showing  their  proper  degeneration  into  simpler 
resultsl  published  earlier. 

Numerical  results  were  computed  and  displayed  for  air-filled  cavities  in 
lossy  rubbers  for  various  cases.  We  initially  considered  cavity  sizes  and  vis¬ 
cous  losses  extending  the  inviscid  results  of  ref.  1.  We  then  used  the  developed 
computer  codes  to  make  additional  runs  for  cases  where  resonating  cavities 
around  (incident)  frequencies  of  20  kHz  and  60  kHz  were  of  concern.  The  analysis 
has  demonstrated  the  dominance  of  the  giant  monopole  resonance  in  all  the  plotted 
cases,  (i.e.,  sound  speed,  attenuations,  real  and  imaginary  parts  of  the  bulk 
modulus,  etc....).  The  cavity  volume  concentrations  were  varied  as  well  as  the 
amount  of  viscous  loss  in  the  absorbing  matrix.  Criteria  for  the  estimation  of 
the  location  of  the  giant  monopole  resonance  emerged  from  the  study  and  we 
observed  in  the  resulting  graphs  how  higher  viscous  levels  shifted  the  resonance 
peaks  and  broadened  their  widths.  Further  analytic  studies  of  these  effects 
should  be  made  in  the  future. 

We  studied  and  discussed  the  static  results  of  Chaban^  and  demonstrated  how 
these  emerge  as  very  particular  cases  of  the  present  results  in  the  absence  of 
viscous  losses,  for  dilute  concentrations,  and  in  the  long  wavelength  limit,  which 
are  the  conditions  for  the  Chaban  results  to  hold  true. 

We  made  a  critical  study  of  the  regions  of  validity  of  the  series  of 
assumptions  behind  the  present  methodology.  We  found  it  universally  valid  under 
the  ordinary  procedures  used  for  high-frequency  absorber  design.  We  were  able  to 
isolate  four  cases  in  which  additional  corrections  are  required  to  modify  the 
present  results.  These  cases  are:  a)  densely  packed  cavities  of  high  volume 
concentration,  b)  large  resonating  cavities  comparable  in  size  to  the  thickness 
of  the  layer  which  contained  them,  c)  shape  corrections  for  other  than  spherical 
cavity  shapes  particularly  when  the  cavities  lie  near  or  on  the  layer  edges,  and 
d)  inclusions  of  solid  (rather  than  fluid)  composition  within  a  lossless  or  even 
a  lossy  matrix.  We  have  developed  a  plan  of  attack  to  analyze  each  of  these 
modifications  in  the  future. 

This  analysis  has  included  a  detailed  study  of  the  effect  of  (non-uniform) 
cavity  size  distributions  on  the  various  effective  quantities  of  interest,  which 
were  quantitatively  displayed  in  many  graphs  (i.e..  Figs.  15-22).  Non-uniform 
(triangular  in  this  case,  c/o  Fig.  14)  cavity  size  distributions  tend  to  have 
the  same  smoothing  effect  on  the  effective  material  parameters  of  the  composite 

i^R.  M.  Christensen,  Mechanics  of  Composite  Materials.  J.  Wilev  &  Sons,  New  York, 
1979,  p.  60. 
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that  seems  to  be  observed  at  higher  viscous  levels.  They  shift  the  resonances, 
reduce  their  amplitudes,  and  broaden  their  widths  over  wider  spectral  regions. 

All  these  considerations  and  quantitative  displays  provide  the  building  blocks  of 
the  trade  off  studies  that  must  be  conducted  by  absorber  designers  in  order  to 
understand  and  construct  efficient  underwater  sound  absorbers. 
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APPENDIX  A  -  PROGRAM  LISTINGS 
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